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INTRODUCTION 


This  interim  report  outlines  the  technical  progress  to  date  ' 
on  the  Electrical  Overstress  Program  at  the  US  Army  Missile  Research 
and  Development  Command  (MIRADCOM) , sponsored  by  Defense  Nuclear  Agency 
(DNA)  , SAMSO  (MNNH) , and  M1RADC(M,  and  provides  a brief  outline  of  the 
background,  plans,  etc.  The  program  Is  aimed  at  providing  a better 
theoretical  model  of  the  physics  of  second  breakdown  In  Junction 
bipolar  diodes,  transistors,  and  Integrated  circuits,  and  applying  this 
Information  In  the  development  of  nondestructive  screening  techniques 
for  eliminating  device  types  and  Individual  samples  which  are  especially 
susceptible  to  high  power  pulses.  Particular  subjects  of  Interest  In 
FY77  Include  short  (submicrosecond)  pulse  effects.  Importance  of  various 
construction  defects,  geometrical  and  structural  effects,  and  applica- 
bility of  existing  Information  to  screening.  The  program  features  a 
continuation  of  DNA  funded  research  using  special  silicon  on  sapphire 
(SOS)  diodes  In  which  the  effects  of  self-heating  are  directly  observa- 
ble, computer  analysis  and  simulation  of  both  conventional  transistor 
and  special  SOS  diode  behavior,  and  selection  and  measurement  of  elec- 
trical parameters  required  for  prediction  of  overstrsss  susceptibility. 

Efforts  currently  under  way  Include  contractor  efforts  In  updating 
and  automating  SOS  diode  test  facilities,  analysis  of  anomalous  short 
pulse  (nanosecond  time  regime)  effects,  generalizing  and  "speeding  up" 
an  existing  computer  coda  for  second  breakdown  simulation,  design  and 
fabrication  of  new  SOS  diode  samples,  development  of  a quantitative 
physical  model  for  filament  formation,  and  application  of  existing  data 
to  an  Improved  engineering-accuracy  prediction  technique  for  failure 
levels  In  conventional  transistors.  Because  these  efforts  have  only 
recently  begun,  results  arc  of  such  a preliminary  nature  that  It  does 
not  appear  advisable  to  report  them  at  this  time.  The  work  discussed 
In  this  report  Is  thus  an  Integration  of  previous  research  (primarily 
by  these  contractors)  and  current  In-house  efforts.  The  in-house  work 
Is  meant  to  provide  approximate  solutions  to  various  problems  which 
enter  Into  analysis  of  second  breakdown,  so  that  these  solutions  can 
serve  as  a starting  point  for  more  sophisticated  numerical  analyses. 
Factors  considered  thus  far  Include  carrier  avalanche  (first  breakdown), 
gsometrlcal  effects,  self-heating,  and  normal  transistor  action;  only 
silicon  devices  are  considered;  and  conditions  under  which  tunneling 
might  occur  are  excluded. 

Second  breakdown  will  be  taken  here  to  mean  the  transition  from 
uniform  to  nonuniform  conduction  In  Junction  bipolar  devices.  It  will 
Include  the  development  of  fllamentatlon  (current  constriction),  but 
not  damage.  ExtresMS  of  operation,  e.g.,  very  short  or  very  long  driving 
pulses,  msy  cause  thermal  effects  to  be  negligible  or  to  dominate,  but 
In  general,  the  process  must  Include  both  thermal  and  electrical 
effects. 


A number  of  phenomena  are  demonstrated  which  fit  the  second  break- 
down criterion  of  producing  a sudden  drop  in  the  voltage  across  a pn 
Junction  which  is  driven  by  a pulsed  constant  current  source.  These 
phenomena  include  nonlinear  heating,  change  in  the  sign  of  the  slope  of 
the  resistivity-temperature  curve,  self-biasing  due  to  lateral  current 
flow  (current  mode),  and  melt  channel  formation.  The  latter  is  excluded 
so  that  the  definition  can  be  further  restricted  to  a nondestructive 
process.  Effects  of  nuclear  radiation  have  not  yet  been  considered 
(the  so-called  synergistic  effects),  but  will  be  Included  in  later 
computer  programming. 

It  has  been  noted,  in  agreement  with  theories  of  a critical  temper- 
ature or  a critical  energy,  that  when  reverse  biasing  current  pulses  of 
less  than  a few  milliseconds  in  duration  are  used  as  the  driving  source, 
the  hot  spots  always  begin  in  the  avalanching  depletion  region.  For 
this  reason,  a large  portion  of  the  work  has  been  directed  toward 
determining  the  physical  location  where  first  breakdown  would  occur. 

A need  for  Information  on  sustaining  voltages  for  nonplanar  structures 
was  clearly  indicated.  Previous  work  by  many  authors  had  already  indi- 
cated that  hot  spots  normally  form  in  avalanche  regions,  but  not 
necessarily  at  the  point  of  maximum  electric  field  as  calculated  from 
the  geometry  of  the  device.  As  shown  schematically  in  Figure  1,  varia- 
tions in  the  geometrical  dimensions  of  a given  device  may  be  sufficient 
to  cause  avalanche  breakdown  to  occur  in  entirely  different  locations 
In  Doelnally  the  same  device  [1].  This  has  important  consequences  for 
second  breakdown,  because  avalanche  is  usually  the  first  step  in  the 
sequence  starting  with  application  of  the  transient  signal,  passing 
through  second  breakdown,  and  ending  in  damage.  One  sample  may  be 
nearly  impossible  to  damage  in  a given  circuit  configuration,  while  a 
second  will  be  extremely  susceptible  to  damage.  It  is  important  to 
note  that  while  second  breakdown,  in  the  reverse-bias  case,  begins  in 
an  avalanching  region,  the  location  of  the  second  breakdown  site  seldom 
coincides  with  the  location  of  maximum  field.  Mlcrostructural  defects 
have  been  proposed  by  some  authors  as  the  explanation  of  this  situation, 
because  they  could  account  for  local  high  field  concentrations.  Other 
authors  disagree.  A large  local  electric  field  could  be  located  at  or 
near  the  hot  spot  if  a mlcrostructural  Inhomogenelty  existed. 

Budensteln,  et  al.  [2]  have  shown  by  calculation  that  the  effect  of  an 
inhostogenelty  is  to  shift  the  hot  spot  location  from  its  geometrically 
predicted  site  to  an  intermediate  location  between  the  site  of  the 
inhomogenelty  and  the  geoisetrlcally  predicted  site.  Thus,  a submlcro- 
scopie  examination  of  the  new  hot  spot  site  would  probably  show  nothing 
out  of  the  ordinary.  Horaover,  the  hot  spot  locations  are  known  to 
vary  with  the  amplitude  of  the  driving  pulse.  Any  attempt  to  correlate 
the  presence  of  an  inhomogenelty  with  an  extemally-maasured  electrical 
parameter  should  be  expected  to  require  driving  the  device  into  avalanche, 
because  the  iahomogenelty  may  very  well  be  dormant  as  far  as  electrical 
activity  Is  concerned  until  the  avalanche  condition  Is  reached.  For 
exaisple,  conditions  required  for  the  generation  of  certain  types  of 
electrical  noise  might  include  driving  the  device  into  heavy  avalanche. 
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to)  AVALANCHE  BHEAKOOWN  FROM  REACH-  (bl  NORMAL  AVALANCHE  BREAKDOWN 
THROUGH  IN  VERTICAL  DIRECTION  FROM  SIDEWALL 

Figure  1.  Cross  section  of  devices  with  nearly  the  same  dimensions. 


SOS  diodes  have  already  proved  to  be  an  extremely  effective  tool 
In  studying  second  breakdown.  It  should  not  be  expected  that  the 
results  of  analyses  of  SOS  diodes  will  carry  over  to  conventional 
transistors,  but  rather  that  the  SOS  diodes  should  act  as  a proving 
groimd  for  the  analytical  tools  developed  for  the  study  of  conventional 
devices.  For  exsaple,  the  heat  flow  patterns  In  the  depletion  region 
of  an  SOS  diode  are  vastly  different  than  those  In  the  depletion  region 
of  an  epitaxial  planar  device,  but  the  same  general  equations  are  used 
In  each  case.  Some  problems  do  exist  with  SOS  diodes.  Surface  and 
Interface  states  are  Important  factors  In  determining  avalanche  break- 
down of  the  samples  [3].  This  complicates  the  analysis  and  leads  to 
conditions  unlike  those  found  In  conventional  transistors.  Other  bulk 
effects  may  be  overshadowed  by  surface  effects.  The  avalanche  voltages 
predicted  by  SOS  diode  analysis  are  greater  than  those  seen  experi- 
mentally. The  extensive  series  resistance  regions  tend  to  mask  some 
purely  Junction  effects  and  overemphasize  bulk  properties.  The  most 
serious  limitation  la  that  except  for  electrical  characteristics,  only 
changes  In  temperature  may  be  observed.  Thus,  Injection  effects  can 
only  be  observed  Indirectly.  However,  It  Is  Important  to  realize  that 
every  major,  and  nearly  every  minor,  observation  of  second  breakdown 
effects  In  SOS  diodes  has  been  found  to  have  a counterpart  In 
conventional  devices.* 

An  analysis  along  the  lines  of  the  Wunsch  model  will  be  given  for 
very  high  amplitude  pulses  of  current  (and  therefore  short  delay  times), 
predicting  the  experimentally-observed  Inverse  relation  between  delay 
time  and  applied  power.  Neither  this  model  nor  the  Wunsch  model  for 
lower  smplitude  pulses,  though  both  are  derived  from  the  heat  equation, 
provides  a physical  explanation  for  the  observed  or  critical  temperature. 
Implicit  In  both  models  Is  the  assumption  that  only  geometrical  effects 
are  Important,  not  statistically-based  mlerostructural  defects.  Exam- 
ination of  experimental  data  shows  thet  the  two  models  predict  general 


*Dr.  P.  P.  Budensteln,  private  cowunlcatlon,  Kerch  1977. 
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trends  around  which  a statistical  data  spread  does  occur.  On  the  other 
hand,  in  experiments  on  devices  in  which  geometrical  tolerances  are 
more  tightly  held,  if  only  by  selecting  all  samples  from  a single  wafer, 
the  statistical  distribution  becomes  quite  narrow.  Thus,  it  can  be 
argued  that  ways  must  be  found  to  elli^nate  the  effects  of  ordinary 
tolerance  variations  between  samples  before  tests  are  run. 

An  analysis  of  the  effects  of  change  in  b'^fc  resistivity  with 
temperature  will  be  given.  This  Indicates  that  the  voltage  drop  across 
the  high  resistivity  side  can  occur  quite  rapidly.  While  this  effect 
is  probably  negligible  in  most  cases  in  conventional  devices,  because 
the  device  must  heat  sufficiently  to  cross  over  the  resistivity- 
temperature  "hump”  while  maintaining  a reasonably  large  undepleted 
collector  region,  it  is  Important  to  the  understanding  of  SOS  device 
phenomena. 

The  melt  analysis  of  Budenstein,  et  al.  [2]  is  repeated,  using  an 
approximation  technique  and  completely  different  boundary  conditions. 

The  two  models  can  be  shown  under  certain  circumstances  to  have  the 
same  general  form,  but  the  relative  Importance  of  heat  loss  and  melting 
terms  can  be  vastly  different.  The  analysis  here  assumes  that  the  melt 
channel  does  not  stop  growing  before  its  diameter  becomes  comparable  to 
the  thickness  of  the  silicon  layer  but  reaches  a ribbon-like  state  and 
then  stops  growing.  Depending  on  actual  physical  data,  the  channel  can 
still  stop  while  cylindrical.  Good  data  on  the  thermal  properties  of 
molten  silicon  do  not  appear  to  exist.  Some  estimates  can  be  made  by 
working  back  from  curve  fits  of  laboratory  data  to  the  theoretical 
equations.  Even  thermal  data  below  the  melting  point  of  silicon  show 
a wide  spread.  This  will  undoubtedly  remain  a problem  when  more 
sophisticated  numerical  analyses  are  performed. 

Finally,  it  seems  appropriate  to  stress  the  preliminary  nature 
of  these  results  and  to  point  out  that  many  of  the  authors  cited  in 
this  report  disagree  strongly  with  each  other's  interpretations  of  the 
data  presented  here. 


II.  TYPES  OF  FIRST  BREAKDOWN 

Avalanche  breakdown  in  pn  Junctions  is  a subject  which  appears 
in  essentially  every  text  on  semiconductor  device  physics,  but  its 
treatment  is  usually  rather  sketchy.  Moreover,  the  solution  for 
cartesian,  or  ractangular,  coordinates  is  the  only  one  discussed,  and 
the  sustaining  voltage,  which  depends  on  a limited  avalanche,  is  seldom 
crasldered;  however,  there  are  the  exceptions  [1,  4,  and  5].  Following 
Gartner  [6],  the  discussion  here  will  begin  with  the  continuity 
equations.  In  each  devlca,  the  volume  of  the  pn  Junction  depletion 
region  can  uaually  be  approximated  by  considering  it  to  be  made  up  of 
subvoltaes  in  the  shape  of  rectangular  parallelepipeds,  quarter  segments 
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of  cylinders,  and  eighth  segments  of  spheres  (Figures  2,  3,  4,  and  5). 
Symmetry  then  allows  these  volumes  to  be  treated  one-dlmenslonally  In 
the  appropriate  coordinate  system,  because  the  current  density  vector 
has  a single  component.  Each  subvolume  considered  will  be  referred  to 
by  Its  coordinate  system:  cartesian,  cylindrical,  spherical.  Because 
the  boundary  conditions  and  Internal  characteristics  (doping  density, 
etc.)  will  not.  In  general,  be  the  same  for  each  subvolume,  first  break- 
down will  occur  In  one  location  at  a lower  applied  voltage  than  In  all 
the  others.  A comparison  of  the  breakdown  of  all  the  subvolumes  thus 
yields  the  first  breakdown.  Analysis  Is  greatly  simplified  If  total 
carrier  currents  rather  than  the  corresponding  current  densities  are 
used. 


The  continuity  equations  for  holes  and  electrons  are 
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Figure  3.  Cartesian  subvolume  of  depletion  region. 


Figure  4.  Cylindrical  subvolume  of  depletion  region. 


Figure  5.  Spherlcel  subvolume  of  depletion  region. 


For  avalanche  conditions,  the  generation  rates  may  be  assumed  to  be 
equal  and  to  be  much  greeter  than  the  recombination  rates: 


r - r ■ 0 
P n 


For  steady  state  conditions.  Equations  (1)  and  (2)  then  become 


7-Jp  - qg 


- -qg  . 


The  standard  assumption  for  generation  rate  in  bipolar  semiconductor 
devices  is  that  a two-carrier  Townsend  avalanche  process  takes  place, 
so  that  the  differential  equation  for  hole  (minority)  current  density  in 
n-type  material  in  cylindrical  coordinates  is  found  to  be 


~ (rJ)"-aJ  - aJ 
r dr  p n n p p 


and  the  differential  equation  for  electron  (majority)  current  density 
under  the  same  conditions  is 


i-;^(rJ)-oJ  +aJ 
r dr  n n n p p 


Integrating  over  the  surface  area. 


— k _ _o  (I  - I ) - o I 
dr  n ' p'  p p 


or,  in  standard  form. 


■dT  - ^“n  - “p^  \ ‘ -“n^  • 

The  ionization  coefficients  in  Equation  (9)  depend  on  electric  field 
intensity  S,  which  must  be  found  by  solving  Poisson's  equation.  Under 
certain  assumptions  to  bo  discussed  later,  the  charge  density  term  in 

Poisson's  equation  la  essentially  independent  of  currents  I and  I , 

P 

so  that  Poisson's  aquation  nay  be  solved  separately.  As  an  example, 
the  npn  transistor  of  Figure  6 is  considered.  For  "ordinary"  avalanche 
breakdown  in  an  npn  transistor,  the  integral  condition  is: 
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In  which  r la  tha  valua  of  r for  which  Equation  (10)  la  balancad.  It 
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will  ba  ahown  In  the  naxc  aactlon  that  under  certain  clrcumataneaa 
(corraapondlng  to  tha  auatalnlng  voltage),  breakdown  can  occur  under 
very  alld  avalanche  condltlona,  with  only  the  requirement  that  [7] 
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Referring  to  Figure  6, 
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In  idilch  r^  la  the  value  of  r for  which  Equation  (12)  la  balanced, 
a c 

Figure  7 la  the  pnp-equlvalent  of  Figure  A.  The  differential  equation 
for  hole  current  la 
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with  boundary  condltlona 
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and  auatalnlng  voltage  condition 
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Numarous  approximations  for  the  variation  of  Ionization  coefficients 
with  electric  field  are  available.  The  most  common  are  [8] 


a • 


-a^/i 


a ■ A^e 


a ■ A-  exp  (a-/&  + a./g  + a.) 


(19) 

(20) 
(21) 


or  combinations  of  these;  the  A and  a constants  are  determined  by  curve 
fit.  Table  1 shows  typical  values  of  A and  a.  Measured  values  are 
shown  In  Figure  8 [8,9]. 


TABLE  1.  COEFFICIENTS  FOR  a CURVE  FIT 
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-1.64688  X 10^ 
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Figure  8.  Meesured  Ionization  coafflclants  In  silicon  [8]. 


Equation  (21)  la  said  to  yield  a maxlaun  of  2Z  error  over  the 
following  ranges  of S [8]: 

2.33  * 10^  < & < 7.47  X 10^  v/ca  for  holes 

1. 69  X 10^  < * < 5.42  X 10^  V/c«  for  electrons 


which  coopara  quits  favorably  to  tha  lOZ  to  20Z  accuracy  quotad  by 
the  ease  rafarence  for  Equations  (19)  and  (20).  The  data  for 
Equation  (21)  In  Rafarence  8,  page  62,  are  reversed;  those  for 
Equation  (20)  are  correct.  Baeausa  a depends  so  strongly  on  &•  the 
accuracy  of  tha  axprasslon  for  a Is  vary  significant  In  calculations 
of  avalancha  breakdown  voltaga;  this  will  ba  dasonstratad. 


Figure  9 ihovs  avalanche  breakdovn  voltage  versus  doping  density, 

calculated  using  the  one-sided  abrupt  approximation  and  Equation  (21) , 
and  Figure  10  shows  the  corresponding  fields. 


Figure  9.  Junction  breakdovn,  one-sided  abrupt,  cylindrical 
and  cartesian  (planar)  coordinates. 


Table  2 shows  a comparison  of  results  using  Equations  (19),  (20), 
and  (21)  for  a typical  case.  Figure  11  shows  the  corresponding  a (r) 

and  ap(r)  using  Equations  (20)  and  (21),  respectively.  The  difference 

a - using  Equation  (21)  has  been  Included  to  demonstrate  that  a 
n p p 

Is  almost  negligible  compared  to  a . The  former  cannot  be  neglected, 

n 

however,  because  It  Is  the  term  which  leads  to  a finite  width  at 
breakdown.  An  Interesting  relation  exists  between  the  ionization 
coefficients  In  ths  range  265  kV/cm  ^ 1 £ 350  kV/cm  (most  of  the  range 
of  Interest  In  avalanche  breakdown) ,~*as  Indicated  In  Figure  12.  The 

approxlmstlon  (\tnlts  of  em'^) 


15 


r 


Figure  10.  Critical  field  for  breakdown,  one-elded  abrupt 
junction,  cylindrical  and  cartesian  coordinates. 


TABLE  2.  SENSITIVITY  OF  RESULTS  TO  ASSUMED  IONIZATION  COEFFICIENT 


14  -3 

(N  • 3 ^ 10  ca  . cylindrical,  one-sided  ^.orupt) 


Equation 

(ua) 

r„  (ua) 
c 

\r 

&(V/ca) 

5 

28.32 

232.6 

3.607  X 10^ 

5 

28.38 

234.1 

3.624  X 10^ 

5 

28.40 

234.6 

3.630  X 10^ 

(22) 

n p 


■ay  be  useful  in  analytical  approslaatlons  for  the  breakdown  voltage. 

A "breakdown"  condition  (unlimited  increaae  in  collector  current 
with  no  increase  in  base  current)  can  occur  in  a co—on  emitter 
circuit  at  collector-base  Junction  voltages  auch  below  the  avalanche 
voltage.  Very  little  avalanching  (M  * 1.02)  is  required.  To  illustrate 
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r Vim) 

Flfur*  11.  lonlMtlon  eo«fflel«&ts,  using  Equations  (20)  and  (21) 


this  condition,  sn  Eb«rs-Mbll  aqulvslant  circuit  of  s pnp  transistor 
with  ths  collsctor-bass  Junction  In  avslanchs  Is  glvan  In  Flgura  13. 
Ths  coUsctor  currant  undsr  avslanchs  Is  glvan  by 


Th«  condition  for  unlimited  increase  in  is 


OfMp-l 


i.e. , such  that  the  denominator  of  Equation  (23)  goes  to  zero. 
The  condition  for  Ig  ■ 0,  Ig  ■ is 


'1  - 0-0  M 

,■  H 

_ 1 - «t"p  J ' 


'n"-CS 


The  sustaining  voltage  is  defined  as  the  voltage  at  which  both  condi- 
tions in  Equation  (2A)  and  (25)  are  satisfied.  The  equation  corraspond- 
ing  to  Equation  (24)  for  an  npn  transistor  is 


aJi  • 1 
t n 


As  noted  previously. 


1-0, 


r 

■I  o e 

I ' 


(o  - o ) dr' 
p n 


1 - . - ; a . 

apa  J 


■ L'--- 


On)  dr- 


in  which  r has  bean  used  as  the  value  of  r-  at  the  austalnlng  voltage 

to  distinguish  it  from  the  value  of  r^  required  for  avalanche.  Figure  14 

shows  the  abrupt,  one-sided  Junction  approslmatlon  for  the 

cylindrical,  spherical,  and  cartesian  geometries.  Figure  IS  shows  the 
effect  of  a different  doping  profile  for  the  cartesian  case  [4]. 
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Figure  IS.  Normalized  sustaining  voltage  versus  current 
gain  for  linear  graded  Junction  transistors  [4]. 


Reachthrough  In  epitaxial  transistors  Is  a special  condition  In 
which  the  collector-base  Junction  depletion  region  becomes  sufficiently 
wide  that  the  entire  lightly-doped  collector  bulk  Is  depleted,  the 
depletion  layer  edge  reaches  the  heavily-doped  substrate,  and  avalanche 
breakdown  occurs  at  an  abnormally  low  voltage.  Expansion  of  the  deple- 
tion region  Into  the  substrate  is  nsgllglbls,  so  a detailed  analysis 
of  the  high-low  Junction  effects  at  the  substrate-collector  bulk  Inter- 
face is  not  necessary.  The  assuogitlon  of  negligible  mobile  carriers  In 
solving  Poisson's  equation  may  still  be  made  for  fields  between  those 
required  for  reachthrough  and  avalanche,  so  that  the  electric  field 
Intensity  may  be  calculated  In  the  usual  maxmer.  Unlike  the  previous 
cases,  the  field  Intensity  at  x^^  (Figure  16)  Increases  from  zero  with 

Increasing  voltage  until  avalanche  breakdown  Is  achieved.  The  field  j 

at  the  Junction  (x  ■ 0)  at  avalanche  due  to  reachthrough  will  be  | 


denoted  The  value  of  Is  found  by  balancing  the  equation 


Figure  16.  Field  variation  In  a substrate-limited  collector. 


A 

X 

/ss  -I  (a  - a ) dx, 

x'  " P 1 

a e 

n 


by  an  Iterative  process,  with  the  upper  limit  fixed.  Similarly,  when 

4(x„)  > 0, 


1 - a. 


(at  - 3 ) dx, 
n p I 


Becauae  the  upper  limit  Is  fixed,  values  of  the  field  at  the  Junction, 
should  be  selected  according  to  the  equation 


qN.x  / qN.x  \ 

& -lAJw  + Y («  , I 

s e ' \ rt  e / 


In  which  values  of  y are  selected  In  the  range  0 ^ y £ 1 and  a curve 
plotted.  Figure  17  shows  * typical  case  [4]. 
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Figure  17.  NomAllred  sustaining  voltage  versus  current  gain 
for  epitaxial  transistors  with  reachthrough  [4]. 


Punchthrough  is  a linking  up  of  the  depletion  layers  of  the 
emitter-base  and  collector-base  junctions.  In  modem  transistors,  the 
base  is  more  heavily  doped  than  the  collector,  so  the  collector-base 
depletion  layer  extends  farther  into  the  collector  than  into  the  base 
for  a given  (reverse)  bias  on  the  collector-base  Junction.  Emphasis 
on  smaller,  faster  devices  has  led  to  very  thin  bases,  thus  Increasing 
the  likelihood  of  punchthrough  and  somewhat  canceling  this  advantage 
of  the  heavily  doped  base.  Becauee  the  emitter-base  junction  depletion 
layer  is  narrower  at  high  emitter  currents  than  it  is  at  low  emitter 
currents,  worst-case  operation  is  likely  to  be  at  high  collector 
voltages  and  low  currents.  Neglecting  any  possible  interaction  effects 
when  the  two  depletion  layers  come  very  close  together,  the  condition 
for  punchthrough  is  simply  found  by  determining  the  width  of  the 

base  side  of  the  collector-base  depletion  layer  (Figure  18)  as  a 
function  of  reverse  bias  voltage  and  comparing  this  width  to  the  distance 
W2  between  the  base  edge  of  the  emitter-base  depletion  layer  and  the 


collector-base  Junction.  Punchthrough  occurs  when  W, 


Avalanche 


i 


Figure  18.  Doping  profile  (schematic). 


multiplication  Is  not  necessary  for  punchthrough,  although  some 
avalanching  may  be  present.  Because  depends  strongly  on  the  details 

of  construction  of  the  emitter  and  base,  greatest  generality  Is  achieved 
In  exhibiting  results  If  Is  kept  as  a parameter.  A practical  design 

or  analysis  problem  would  then  be,  "Given  an  epitaxial  double-diffused 
transistor  with  a punchthrough  width  W^,  a collector-base  junction  side- 


wall radius  rj,  and  a collector  epitaxial  layer  thickness  d,  which  would 

occur  first,  avalanche  (either  sidewall  or  vertical),  reachthrough,  or 
punchthrough?"  Fortunately,  Is  already  calculated  as  x^^  - Xj  In 


finding  the  relationship  between  the  width  of  the  base  side  of  the 
collector-base  depletion  layer  and  that  of  the  collector  side.  Results 
of  such  a calculation  are  shown  In  Equation  (21)  and  Figure  15  for  a 
Gaussian  diffusion  Into  a constant  epitaxial  doping,  roughly  compatible 
with  a Hotorola  2N2222.  Figure  19  shows  reachthrough  width  versus 
collector-base  bias  voltage  [8].  Phillips  [10]  points  out  that 


(32) 


In  which  x^  Is  the  width  of  the  depletion  region  at  ordinary  breakdown. 
This  can  be  rewritten 
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Figure  19.  Avalanche  and  reachthrough  voltages  [8]. 
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rt 


= d 


'2qB  V, 


(33) 


For  a Motorola  2N2222,  epitaxial  doping  density  ■ 10^^  cm  ^ , and  a 

typical  value  of  d Is  8 ym.  Equation  (33)  predicts  • 202,  In  good 

agreement  with  Figure  19.  It  should  be  noted  that  If  d - 8 um,  then 
■ 0.75  um.  Because  the  base  width  Is  nominally  1 um,  reachthrough 

will  probably  occur  before  punchthrough  at  normal  forward  bias. 


III.  SECOND  BREAKDOWN 

The  literature  on  second  breakdown  Is  sufficiently  volisslnous 
that  no  attempt  will  be  made  to  susmarlze  more  than  a few  key  papers 
here.  Budenstaln  [11]  made  a survey  of  papers  published  before  1970, 
Schafft  and  French  [12]  made  a survey  of  very  early  work,  and  the  IEEE 
Transactions  on  Electron  Devices  devoted  two  special  Issues  to  the 
subject  [12,13].  These  three  references  can  be  used  as  a source  for 
nearly  all  of  the  early  work  on  second  breakdown.  The  moat  recent 
published  work  In  the  area  at  the  time  of  this  writing  Is  probably 
that  of  Hard  [14],  Sutherland  and  Kannady  [15],  and  Rabum  and  Causey 
[16,17]. 
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Early  work  Involved  Che  search  for  a crlCical  parameter  upon  which 
CO  base  analyses  of  second  breakdown.  Ford  [18]  and  Melchior  and 
Struct  [19]  suggested  a critical  temperature,  defined  to  be  that  temper- 
ature at  which  Che  Intrinsic  carrier  concentration  reached  the  doping 
density  on  the  lightly-doped  side  of  the  pn  Junction.  Schafft  and 
French  [20]  proposed  a critical  energy,  based  on  Che  fact  that  a delay 
time  occurs  before  breakdown  when  a pulse  of  constant  power  Is  applied. 
Delay  time  to  breakdown  should  thus  be  Inversely  proportional  to  Input 
power;  this  Is  seen  at  very  high  power  levels  for  which  delay  times  are 
less  chan  100  nsec.  Wunsch  and  Bell  [21]  used  the  critical  temperature 
model  to  derive  an  expression  predicting  that  the  square  root  of  the 
delay  time  would  vary  Inversely  with  Input  power.  This  model  bridges 
the  time  gap  between  the  simple  Inverse  dependency  observed  only  below 
a few  nanoseconds,  and  dc  breakdown.  They  assumed  a critical  tempera- 
ture T^  for  Junction  failure  and  Idealized  the  Junction  heat  source  as 

a plane  source  of  area  A (which  might  be  a constricted  region) , located 
at  the  origin  of  coordinates.  The  previous  chapter  Illustrates  the 
complexities  of  even  a slightly  more  detailed  model.  The  following  Is 
a simplified  explanation  of  their  derivation.  The  heat  equation  for 
Joule  heating  and  negligible  thermal  diffusion  Is 


P 


m 


(34) 


In  which  p is  Che  mass  density  and  c Is  Che  heat  capacity  or  specific 
m 

heat  per  unit  mass.  Integrating  over  Che  volume, 

p c FA  - P (35) 

n 

which  defines  a heat  source  F In  terms  of  Che  dissipated  power  P and 
area  A: 


F 


P 

A 

IB 


(36) 


The  source  ftincclon  for  a plane  source  at  the  origin  is 


(37) 
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so  that,  convolving  F and  G, 


P/tJ 


P_  A 

ID 


^TTQI 


(38) 


In  which  t^  Is  the  time  between  application  of  an  electrical  pulse  and 

noticeable  change  in  electrical  characteristics;  e.g.,  application  of 
a large  "rectangular"  pulse  of  constant  current  leads  to  a roughly 
constant  voltage  until  time  t^,  when  the  voltage  suddenly  drops. 

Rearranging  Equation  (38),  the  "Wunsch  model"  Is  found: 


P 

A 


c /iro  (T-  - T ) 
m t * 


const 


(39) 


The  model  predicts  the  general  trend  in  Junction  behavior  over  a large 
-8  -2 

time  regime,  10  S t^  S 10  sec  for  most  devices,  but  (as  noted  by 

Wunsch  and  Bell)  does  not  explain  the  large  variation  about  the  predicted 
value  of  P for  any  given  t^.  The  problem  of  sidewall  avalanche  break- 
down versus  vertical  reachthrough  discussed  In  Section  II  Is  one  step 
In  explaining  this  variation. 

The  critical  energy  model  may  be  used  to  demonstrate  the  effects 
of  short,  high  current  pulses  during  which  the  heat  generated  has  no 
time  to  leave  the  depletion  layer.  Whether  current  constricts  during 
the  rise  time  of  the  pulse  or  not  at  all  Is  actually  unknown.  An 
Idealized  analysis  will  now  be  given  which,  like  the  Wunsch  model.  Is 
very  approximate  but  gives  some  Insight  Into  the  effects  of  various 
parameters.  It  Is  assumed  that  the  current  constricts  to  a channel  or 
channels  of  constant  total  area  A^  (which  will  be  some  fraction  of 

Junction  area  A)  and  length  L In  a ona-sldad  abrupt  Junction  of  doping 
density  N.  Because  A^  will  actually  vary  with  time,  the  analysis  will 

only  be  valid  for  pulses  short  enough  that  A^  does  not  vary  significantly. 

The  value  of  L can  be  approximated  by  the  distance  from  metal  contact 
to  Junction  on  the  lightly  doped  side,  or  photographic  data  can  be 
used.  It  Is  assumed  all  current  flows  through  A^.  Then 

J - ^ (40) 
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and 


o L ^2 

The  heat  generated  per  unit  volume  within  the  constrlctlon(s)  may  be 
expressed  as 


2 I V, 


J-&  - 


The  subscript  J has  been  used  on  voltage  and  power  to  emphasize  that 
these  are  values  for  the  junction  region,  rather  than  external  values. 
The  heat  equation,  assuming  an  average  value  for  K,  luy  then  be  written 

^2,  2 

3T  9 T , J / A ^ \ 

3?  ® “T  2 • 

3x^  c p A.L 

ID  t 

The  boundary  conditions  which  will  be  assumed  are 


T(x,0)  - T 


T(0,t) 


U (L.t)  - 0 . (A6) 

Equation  (45)  is  probably  not  too  bad  an  assumption  for  the  short  pulses 
of  Interest.  The  last  boundary  condition  [Equation  (46)]  Is  obtained 
by  considering  a synnetrlc  abrupt  junction  for  which  the  peak  tempera- 
ture Is  assuiDSd  for  simplicity  to  occur  at  the  metallurgical  junction. 

A numerical  solution  to  this  problem  can  be  found  by  using  standard 
techniques;  however,  due  to  the  approximations  made,  a simpler  solution 
technique  seems  warranted.  If  e solution  of  the  form 


n 

T • T^  ^ aj^(t)x^ 


Is  assumed  using  boundary  conditions  In  Equations  (14) , (15)  and  (16) 
and  n ■ 2,  then 
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For  small 


o P„  1.  - I.) 
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Pa  - T.) 
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Using  n “ 3, 


T,  - T(L,t-)  - T + 


“ 3c  p A,  a 
n r 


(i . 


For  small  t^.  Equation  (49)  is  again  obtained.  Equation  (49)  may  also 

be  obtained  by  using  Laplace  transforms.  Solving  Equations  (44)  through 
(46)  leads  to  a complicated  expression  which  may  be  simplified  by  the 
fact  that  the  solution  for  small  time  t corresponds  to  large  s In  the 
transformed  solution.  This  approximation  again  leads  exactly  to 
Equation  (49) . Equation  (49)  may  be  rewritten  in  terms  of  externally- 
measured  values  and  voltage  drop  across  the  ohmic  bulk 

region,  as  ^ 


* 


, - V.  ,,  ) A- 
>1  bulk  f 


in  which  P is  the  total  power  dissipated,  and  is  the  applied 

voltage  (assumed  constant).  If  negligible,  a simple  Inverse 

relation  exists  between  P and  t^.  If  not,  an  apparent  "data  spread" 

will  occur.  The  same  argument  obviously  applies  to  the  Wunsch  model. 
Some  other  points  should  be  made  about  Equation  (49).  A single,  con- 
stant T^  is  assumed,  as  in  the  Wunsch  model;  only  the  qualitative 

success  of  the  Wunsch  model  supports  such  an  assumption.  Experimental 
evidence  exists  for  an  inverse  relationship  between  and  t^,  but 

experimental  difficulties  are  encountered.  Equation  (51)  predicts  a 
dependence  on  doping  density  N.  The  effect  of  the  square  root  Is  to 
lessen  the  data  spread  due  to  doping  variations  between  samples. 

A technique  for  direct  observation  of  temperature  variation  In 
second  breakdown  was  developed  by  Sunshine  [22]  and  later  used  by 
Oumln  [23].  This  technique  gave  visual  proof  of  the  current  constric- 
tion theory  proposed  by  many  authors  [12,13].  Budensteln  and  his 
students  [2,24]  later  refined  this  technique  and  made  an  extensive 


study  of  current  constriction  or  fllamentatlon  which  forms  the  back- 
ground for  much  of  this  report.  The  following  is  a brief  description 
of  some  of  Budenstein's  work. 

SOS  diode  test  structures  (Figure  20)  can  be  made  sufficiently  thin 
that  light  can  pass  through  both  the  silicon  layer  and  the  sapphire 
substrate.  The  light  transmissivity  at  any  point  in  the  structure  is 
proportional  to  the  local  temperature  of  the  silicon.  Using  a lateral 
Junction  which  is  viewed  edge-on  (l.e.,  perpendicular  to  current  flow) 
has  a number  of  advantages:  Current  constrictions  are  easily  observable 
(as  darkened,  high  temperature  areas),  reachthrough  problems  of  a 
vertical  Junction  are  avoided,  and  silicon-sapphire  Interface  effects 
are  less  significant  [25].  The  limitation  of  the  technique  is  that 
only  the  effects  of  self-heating  are  directly  observable.  Figures  21 
and  22  show  the  optical  arrangement  used  by  Budenstein,  et  al.  [2,24] 
to  make  temperature  measurements.  The  probes  shown  in  Figure  22  are 
used  to  pulse  the  device  and  to  measure  voltage  and  current.  Figure  23 
shows  photographic  results  using  this  arrangement.  Constant  current 
pulses  were  used,  and  sufficient  time  was  allowed  to  pass  for  develop- 
ment of  well-defined  hot  spots  (the  dark  spots  on  the  photograph)  at 
the  Junction.  A large  (reverse-bias)  current  pulse  leads  to  the  forma- 
tion of  many  hot  spots,  while  lower  currents  lead  to  fewer  hot  spots. 

The  uniformity  of  spacing  is  the  rule,  rather  than  the  exception. 
Sonuniformlty  is  probably  related  to  local  inhomogeneity.  As  shown 
in  Figure  24,  the  saturation  (maximum)  number  of  hot  spots  in  the 
particular  device  used  increased  roughly  linearly  with  the  amplitude 
of  the  constant  current  pulse  for  currents  greater  than  approximately 
200  mA.  If  each  hot  spot  is  sufficiently  isolated  from  the  others 
that  it  can  be  treated  Independently  and  if  a critical  current  density 
exists  for  onset  of  fllamentatlon,  then  the  experimentally  observed 
uniformity  of  hot  spot  size  [2,24]  suggests  that  the  number  of  hot  spots 
should  Increase  linearly  with  applied  current.  "If  pulse  duration  is 
kept  constant  while  progressively  Increasing  pulse  amplitude,  the 
number  of  filaments  increases  to  a saturation  value.  Further  increase 
in  pulse  amplitude  simply  drives  more  current  through  existing  filaments 
and  ultimately  leads  to  melt  formation"  [2].  Very  large  reverse  bias 
pulsing  leads  to  narrower  filaments  (elongated  hot  spots),  and,  as 
predicted  by  the  model  just  discussed,  more  of  them  than  at  lower 
current  levels.  Forward  bias  pulsing  leads  to  a single  hot  spot  which 
is  much  wider  (more  diffuse)  than  those  seen  at  reverse  bias  pulsing. 

The  effective  cross  section  of  this  hot  spot  is  greater  for  current 
conduction  and  heat  sinking  then  that  of  a single  hot  spot  caused  by 
reverse  bias  pulsing,  so  It  is  not  surprising  that  damage  occurs  at 
higher  currents  and  lower  voltages  than  those  required  for  damage  at 
reverse  bias.  An  extensive  analysis  of  forward  bias  second  breakdown 
has  been  performed  by  Navon  [26]. 

In  the  reverse  bias  case,  filament  formation,  or  hot  spot  elonga- 
tion, proceeds  as  illustrated  schematically  In  Figure  25.  The  width 
changes  very  little  as  the  length  Increases.  The  extemally-measured 
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Figure  20.  SOS  diode  geometry.  Doping  materials  are  diffused 
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Figure  22.  Experimental  arrangement  of  the  optical  system 
used  In  making  stroboscopic  observations  [24]. 


voltage  should  decrease  more-or-less  linearly  with  the  elongation  of 
the  filament  and  may  account  for  the  gradual  decrease  In  voltage  seen 
In  some  voltage  waveforms.  When  the  filament  reaches  the  other  side 
of  the  high  resistivity  region,  the  "ballasting"  effect  of  the  high 
resistivity  material  is  gone  and  essentially  all  of  the  voltage  drop  is 
across  the  filament,  which  is  already  so  hot  that  the  material  In  the 
filament  is  at  or  near  the  intrinsic  temperature.  A melt  channel  forms 
within  the  intrinsic  channel  almost  immediately  after  the  high  resis- 
tivity region  is  bridged.  An  interesting  Indication  of  the  temperature 
distribution  around  the  hot  spots  is  given  schematically  in  Figure  26 
which  shows  the  decrease  in  avalanching  next  to  the  filament,  due  to 
lateral  heat  flow.  A region  of  normal  avalanche  is  seen  in  close 
proximity.  The  sapphire  substrate  acts  as  a heat  sink  not  found  in 
conventional  diodes  and  transistors,  so  the  latter  are  expected  to  have 
I somewhat  different  thermal  profiles  than  those  shown  in  Figure  26. 

I The  filament  on  the  right  in  Figure  26  is  growing  at  the  expense  of 

i the  others:  it  will  soon  rob  most  of  the  current,  stunting  the  growth 

rate  of  the  other  filaments  and  increasing  its  own.  Figures  25  and  26 
j are  based  on  observations  made  while  watching  a time-lapse  color  movie 

film  of  second  breakdown  in  SOS  diodes.  The  film  was  produced  by 
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Figure  24.  Saturation  number  of  hot  spots  versus 
critical  current  for  breakdown. 


Budensteln,  et  al.  [2]  and  has  been  shown  at  several  technical  meetings. 
Figure  27  shows  a thermal  mapping  along  a line  (Indicated  by  horizontal 
arrows)  parallel  to  and  at  a separation  of  12.7  pm  from  a junction  In 
which  the  hot  spots  are  clustered,  rather  than  uniformly  spaced  [24]. 

Two  very  small,  v-shaped  metal  protuberances  (Indicated  by  vertical 
arrows)  along  the  edge  of  the  upper  metal  pad  seem  to  have  a significant 
effect  on  the  location  of  the  hot  spots.  The  protuberances  and  the 
surrounding  high  temperature  darkened  regions  are  far  more  evident  on 
the  original  photograph.  The  protuberances  were  not  seen  Initially  and 
were  observed  to  grow  during  the  multiple  pulsing  required  In  setting 
up  and  taking  data,  possibly  by  electromigration  of  aluminum  from  the 
metal  pad.*  Because  the  location  of  the  hot  spots  did  not  change  over 
a long  period  of  time  (Each  dot  on  the  graph  In  Figure  27  represents  a 
statistical  average  of  hundreds  of  pulses,  so  that  a week  was  required 
to  obtain  the  data  for  this  single  graph.).  It  may  be  that  nearly  all 
of  the  electromigration  took  place  very  quickly.  This  would  be  con- 
sistent with  a diffusion  process.  It  should  be  mentioned  that  the 
averaging  process  Is  known  to  have  lowered  the  "peaks"  and  filled  In 


*?.  P.  Budensteln,  private  communication,  November  1976. 
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Figure  26.  Temperature  distribution  along  junction  (schematic). 

the  "valleys"  in  the  thermal  map  of  Figure  27.*  Such  an  effect  of 
repeated  pulsing  Is  of  more  than  academic  Interest;  overvoltage 
transients  are  often  In  the  form  of  damped  sine  waves,  so  the  device 
could  be  hit  several  times  In  rapid  succession. 

The  work  of  Berglund  and  Klein  [27]  serves  to  round  out  the 
Information  on  combined  thermal  and  electronic  effects  by  providing 
a survey  of  thermal  effects  on  switching  of  solids  In  general  li. 
points  out  that  the  discharge  of  the  Junction  capacitance  onxsc  be 
conaldercd  In  a thorough  transient  analysis  of  second  breakdown.  This 
discharge  current  is  the  displacement  current  term  in  the  model  used 
by  Rabum  and  Causey  [16]. 


1 


I 


r 

^ In  SOS  diodes  and  In  some  conventional  transistor  structures,  the 

^ change  In  resistance  with  time  before  current  constriction  may  be 

significant.  A uniformly  conducting  SOS  diode  of  cross-sectional  area 
I A Is  considered.  Because  the  standard  test  condition  Is  a pulse  of 

constant  current  1, 


V(t)  - IR(t)  - p(t)  . (52) 

Considering  only  heat  Loss  to  the  substrate  and  neglecting  heat  flow 
from  the  depletion  region  Into  the  bulk,  such  that  effectively  only  a 
resistive  sheet  Is  considered  and  diffussion  is  negligible,  the  heat 
equation  becomes 


CP  ^ - J-i  - H(T  - T ) - p - H(T  - T ) . (53) 

m dt  * .2  » 

A 

The  resistivity  Is  given  by  the  approximation 
b T 

pj^  e T < 523*K  (54) 

p - p.  - a(T  - T )^  523®K  < T < 623*K  (55) 


[ P2  e ^ 623*K  < T (56) 

In  which  Pq,  Pj^,  P2,  a,  b,  and  c are  found  by  curve  fitting,  and  Is 

the  temperature  at  the  peak  of  the  resistivity  curve.  The  Initial  tem- 
perature Is  assumed  to  be  T^  • 300*K.  Equation  (53)  may  be  solved  as 

a nonlinear  first  order  differential  equation  by  using  a Runge-Kutta 
technique  and  the  appropriate  equation  for  the  resistivity  in  each 
temperature  regime.  Since  the  right-hand  side  of  Equation  (53)  is  a 
function  only  of  T,  direct  integration  is  possible,  so  that  for 
T < 523*K, 


for  523"K  < T < 623*K 


(/?"  ^)(/t - ^ * 0 -(/t  - ^)(/?^  - ^)« 


ce_  A 


2^  1“ 

—V  - *1> 

ce_  A 


(c  - tj 


(58) 


in  which  T(tj^) 


Tj^  - 523*K,  and  for  623*K  £ T, 


b2T2  . ^ ‘^2^2 

® 12 

CP  A 

tD 


(59) 


in  which  T(t2)  ■ T2  ■ 623*K.  The  corresponding  voltage  in  each  case 


is  found  by  calculating  the  resistivity  and  substituting  it  into 
Equation  (52).  For  example.  Equation  (57)  would  be  substituted  into 
Equation  (54),  and  Equation  (54)  would  be  substituted  into  Equation  (52). 
Substituting  some  typical  data  [11],  the  effect  of  self-heating  over  the 
resistivity  peak  while  still  in  uniform  conduction  is  found  to  be  as 
shown  in  Figures  28,  29,  and  30.  This  may  explain  some  of  the  pertur- 
bations seen  on  the  voltage  waveforms  [11]  prior  to  the  sudden  transi- 
tion to  a lower  voltage  state.  Sutherland  and  Kennedy  [15],  Rabum  and 
Causey  [16],  and  Budensteln,  et  al.  [2]  have  used  a thermal  equivalent 
circuit.  Their  results  Indicate  that  such  an  analysis  is  feasible  if 
the  source  of  the  initial  Inhomogenelty  can  be  found.  As  an  example 
of  such  a thermal  equivalent  circuit.  Figure  31  shows  the  Sutherland 
and  Kennedy  model  [15].  represents  the  power  dissipation  driving 

term  in  the  heat  equation,  passive  components  R^,  R.^,,  and  represent 


terms  arising  from  finite  difference  approximation  to  the  partial 
derivatives,  and  represents  the  substrate  temperature  (assumed 

invariant,  l.e.,  an  ideal  heat  sink). 


IV.  MELT  CHANNEL  FORMATION  AFTER  SECOND  BREAKDOWN 

The  conditions  leading  to  second  breakdown  in  a pn  junction 
can  be  very  difficult  to  analyze,  even  with  a digital  computer,  but 
once  current  constriction  or  plnch-ln  has  occurred,  the  analysis  of 
actual  damage  can  be  greatly  simplified.  Budenstein,  et  al.  [2],  have 
shown  by  direct  observation  that  damage  due  to  second  breakdown  in  a 
pn  Junction  is  the  result  of  a melt  channel  forming  within  the  intrinsic 
("wash-out”)  pinch-in  channel  across  the  depletion  region.  Basically, 
melt  channel  formation  is  a moving-boundary  problem  with  phase  change 
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Figure  30.  Temperature  variation  around  resistivity  peak. 


and  cylindrical  syrmtatry.  Computer  solutions  for  this  problem  under 
a variety  of  assumptions  [2,25,28]  have  existed  for  some  time,  but  an 
approximation  technique  of  easy  applicability  has  not.  Moreover,  some 
questions  exist  as  to  the  assumptions  made  In  these  and  other  earlier 
analyses.  The  analysis  of  a related  planar  problem  by  Goodman  [29] 
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will  be  compared  (but  Finlayson's  notation  [30]  will  be  used),  particu- 
larly with  reference  to  boundary  conditions.  In  the  simplistic  model 
to  be  used  here  (Figure  32),  a single  melt  channel  Is  assumed  and  Is 
represented  by  a cylinder  of  radius  s(t)  and  length  L,  surrounded  by 
a large,  ellipsoidal  Intrinsic  region  or  chatmel.  The  total  current 
I will  be  a constant  (Figure  33),  following  Budensteln's  experiments 
and  Is  given  by 


J(r,t)2irrdr 


v^ti 

L 


2iirdr 


V(t)  - 


a(r)2Trrdr 


In  which  J(r)  Is  flowing  In  the  z direction  but  varies  In  Intensity 
with  radial  position.  Current  flow  outside  the  melt  channel  will  be 
neglected  for  the  moment  but  will  be  considered  later.  The  object  of 
this  analysis  Is  to  calculata  the  latter  portion  of  the  V(t)  curve, 
i.e.,  after  the  melt  channel  begins  to  form.  The  first  part  of  the 
solution  will  not  be  accurate,  because  the  assumption  of  negligible 
flow  outside  the  melt  channel  Is  poorest  then.  For  definiteness,  the 
configuration  shown  In  Figure  34  will  be  analyzed  and  results  compared 
with  those  of  Budensteln,  et  al.  [2],  page  128.  The  general  equation 
Is 

= ».  If  • s <■= ».)  It]  • 

If  K,  c,  and  p are  assumed  to  be  Independent  of  T (e.g..  If  average 
s 

values  are  used),  the  term  In  square  brackets  disappears.  If  this 
assumption  Is  not  made,  it  Is  still  possibla  to  let 


•4/ 


In  which  T is  the  melting  tasverature  of  silicon,  and  K Is  the  value 

■ m 

of  K just  shove  the  melting  tampersture.  Then 


and  tha  Fourlar  aquation  ac  r • 0 


) 


11 

3r 


r"0 


0 


(70) 


(which  la  raqulrad  for  a continuous  sacond  darlvatlva) . Tha  symbol 

danotas  tha  solid  slda  of  tha  malt'^solld  intarface.  Latant  haat  of 
malting  X in  J/kg  antars  only  at  tha  boundary:  tha  rast  of  the  melt 
channel  Is  above  tha  malting  tasiparature. 

If  a surface  (a.g.,  tha  substrata  or  tha  top  of  the  chip)  loss 
term  Is  Included,  Equation  (67)  can  be  written 


e p 


II 

m 3t 


H (T  - T,)  + g 


(71) 


Quantities  o,  p^,  K,  X,  and  c are  taken  to  be  constants  and  equal  to 

thalr  astlmatad  mean  values  over  tha  temperature  range  of  Interest 
(Table  3).  Where  a range  of  temperatures  Is  Indicated,  the  numbers  In 
parentheses  were  used  for  the  calculations.  Their  approximation  will 
be  justified  later.  Than  Equation  (64)  can  be  written 


and  Equation  (64)  can  be  written 


IL 

2 

ITS  a 


In  which 


(73) 


K 


cp_L 


1^0 


2 4 

CP  V s o 


it^oKs^ 


const 

4 


(74) 

(75) 


For  simplicity,  tha  difference  between  R above  and  below  the  melting 
point  Is  neglected. 


■ 
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TABLE  3.  ESTIMATED  THERMAL  PROPERTIES  OF  MOLTEN  SILICON 


Symbol 

Parameter 

Measurement* 

**m 

Mass  Density 

2490  kg/m^ 

c 

Heat  Capacity  per  Unit  Macs 

917  to  1050  (982)  J/kg-*K 

X 

Latent  Heat  of  Malting 

9.17  X 10^  to  1.8  10® 

(1.8  X 10®)  J/kg 

K 

Thermal  Conductivity 

22  W/m-*K 

o 

Electrical  Conductivity 

10®  l/ohm-m 

a 

Thermal  Dlffuslvlty 

9 X 10  ® m^/sec 

T 

m 

Melting  Temperature 

1683  to  1696  (1683)  *K 

^b 

Bolling  Temperature 

2628  to  3418*K 

*Th«  tiser  Is  urged  to  be  cautious  of  these  figures,  which 
are  taken  from  a variety  of  sotirces. 
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dt 


(76) 


r»s 


-0-f 

3r 


r*s 


4*  + il 

dt  3t 


r“s 


3r 


Similarly,  on  the  other  side  of  the  melt-solid  Interface, 


- 0 - il 

-^1 

(r  ' 

° 3? 

dt  r 

3r  ' 

V 3r/ 

• 

!•+ 

The  two  expressions  for  the  Interface  velocity  must  be  equal: 


(77) 


dt  dt 


(78) 


Solutions  to  the  heat  equation  must  be  found  for  each  side  of  the 
Interface.  If  an  approximate  solution  of  the  form 


T(r,t)  • T,  - «(•  - r)  + b(s  - r)^ 


(79) 
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ia  aasumad  for  tha  male  alda,  than 


da  ^ 2a  g 

dt  a a 


(80) 


To  find  an  axpraaalon  for  tha  tamparatura  gradlant  outalda  the  malt, 
it  la  nacaaaary  to  aolva 


il  . ® _L  ^ 3T\ 
3t  " r 3r  3r  / 


(81) 


aubjact  to  tha  condltlona 


T(r,0)  - 

(r  > s) 

(82) 

T(»,t)  - 

(83) 

T(s,t)  - 

(84) 

In  which  la  tha  temperature  at  the  edge  of  the  hot  apot.  The  hot 

spot  la  so  large  compared  to  the  melt  channel  chat,  for  the  duration 
of  a typical  test  pulse.  It  may  be  considered  Infinite  In  size  and 
Initially  at  constant  temperature  In  the  raglon  around  the  malt.  This 
la  Juaclflad  by  tha  fact  that  the  tamparatura  peaks  In  Figure  25  are 
on  the  order  of  SO  ym  In  width,  while  tha  solution  hare  la  Intended  to 
be  valid  only  until  the  malt  diameter  reaches  approximately  1 ym.  A 
more  sophlselcatad  approach  to  the  SOS  diode  analysis  would  have  to 
Involve  the  assumption  of  an  Initial  temperature  distribution  Inside 
Che  hoc  spot  (which  would  also  have  an  outer  moving  boundary)  and  the 
generation  of  heat  Inside  the  Intrinsic  channel  due  to  the  small,  but 
finite,  current  flowing  there.  For  a conventional  transistor,  o^y 
Che  substrate  loss  term  could  be  neglected.  Continuing  with  the 
analysis. 


T - + (T^  - T^) 
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(86) 


in  which  Ej^(y)  - -El(-y)  is  ths  Incsgral 

m 

Ei(y)  -y*  ^ dt  . 
y 

From  Equation  (77), 


(87) 


Tha  tasparatura  (r  £ s)  may  Chan  ba  wrlccan 


T ■ T -f  His — ? 
n 2 

8 + 4at 


(88) 


Equation  (85)  may  ba  usad  to  find  T for  r > s,  if  desired.  Equations 
5.1.53  and  5.1.56  of  Abramowltz  and  Stegun  [31]  may  be  used  to  evaluate 
the  integral  in  Equation  (86).  Tha  expression  for  s(c)  is  found  from 
either  Equation  (80)  or  Equation  (87) : 


8 ■ 8. 


(89) 


Equation  (88)  chan  becomes 


T • T + 


1 (s^  - r^)  . ^ . 

2 2 A (if  i ®) 

Lt^oIC(«^  + 4a)  J s 


and  Equation  (85)  bacomas 


T - T.  + (T_  - T.) 


ilS) 

■.«) 


(r  i a) 


For  the  data  of  iotarast  hare. 


2in  + Y - 

T.  -f  (I  - T.)  * L ■?■»•/  fl 

* J M \ \2 


r > a 


(90) 


(91) 


(92) 


4« 


in  which  y - 0.57721... 
voltage  Is 


is  Euler's  constant.  From  Equation  (67),  the 
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Figure  35  shows  the  melt  channel  radius  versus  time  for  an  SOS  diode. 

The  physical  model  breaks  down  at  very  short  times  such  that  the 
resistance  of  the  melt  channel  Is  comparable  to  the  bulk  resistance. 

As  the  melt  channel  diameter  becomes  comparable  to  the  thickness  of 
the  silicon  layer.  Its  cross  section  will  no  longer  be  circular  and 
the  model  Is  again  Inaccurate.  Similar  arguments  hold  for  an  ordinary 
diode  or  transistor  Junction.  In  a small  transistor,  the  channel 
diameter  could  become  comparable  to  the  channel  length,  so  that  neglect- 
ing heat  flow  from  the  ends  of  the  channel  would  no  longer  be  valid. 

For  the  SOS  diode  shown  previously  In  Figure  34,  and  I ■ 0.56  A,  the 
temperature  la  essentially  Independent  of  position  throughout  the  melt 
region.  In  other  words,  the  second  term  In  Equation  (90)  contributes 
almost  nothing.  The  derivative  of  T with  respect  to  r Just  Inside  the 
Interface  Is  not  negligible  and  dominates  the  growth  rate  of  s(t). 

Figure  36  shows  the  temperature  variation  In  and  around  the  melt  channel 
In  an  SOS  diode.  The  particular  SOS  diode  used  for  this  example  was 
relatively  heavily  doped,  so  that  the  partially-avalanching  depletion 
region  extended  only  a small  portion  of  the  way  along  the  melt  channel. 
The  Intrinsic  channel  was  therefore  surrounded  by  and  In  series  with 
ohmic  bulk  material.  The  temperature  profile  corresponds  to  t^  of 

Figure  37,  which  shows  profiles  for  an  epitaxial  transistor,  e.g.,  a 
2M2222.  Voltage  variation  with  time  depends  strongly  on  resistance 

of  the  material  In  parallel  with  the  melt  channel  and  will  be  considered 
later. 

A very  desirable  laq>rovement  can  be  made  by  considering  the  current 
outside  the  melt  channel  to  flow  through  a constant  reslstanca  (this 

Is  a good  approximation  [2];  nearly  all  of  this  current  will  flow  through 
the  Intrinsic  channel.  The  heat  generation  by  this  external  current 
can  bo  shown  to  be  negligible  compared  to  that  generated  Inside  the 
melt.  Then  the  melt  channel  and  the  rest  of  the  Junction  share  total 
current  I (Figure  38)  and 
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Figure  37.  Temperature  profile  around  melt, 
collector  of  epitaxial  tranalator. 


in  which  ia  the  bulk  reeiatance  in  parallel  with  the  melt  channel. 

The  radiua  of  the  intrlnalc  channel  increaaea  with  that  of  the  melt 
chaxmel,  because  the  region  around  the  melt  channel  becomes  hotter  as 
the  melt  channel  grows.  At  t * t^. 


Vo  - V(t,) 


1 ■ ^^*0 
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r 


V(t) 


Figure  38.  Equivalent  circuit. 

In  which  8(t  ) ■ S-. 

o 0 

For  Rg  large  coopered  to  the  melt  channel  resistance,  which  will 

be  true  for  all  but  very  short  times,  and  s(t^)  very  small,  the  solution 

obviously  approaches  that  of  the  simpler  model.  The  function  g(s)  now 
becomes  far  more  complicated 


(96) 


From  Equations  (78),  (80),  and  (87),  It  may  be  seen  that  the  temperature 
gradient  at  r • s (or,  equivalently,  the  parameter  a)  adjusts  with  heat 
generation  such  that  s(t)  Is  Independent  of  R^.  However,  V and  T depend 

on  through  g.  As  noted  previously.  Equation  (88)  for  the  temperature 

In  the  melt  Is  only  significant  In  that  It  can  be  used  to  calculate  the 
te8q>eracure  gradient,  which  depends  linearly  on  g.  Using  Equation  (96), 
the  temperature  gradient  at  s for  finite  external  resistance  R^  Is 
given  by 

3T|  2s»  -2sq 

+ «o  ■ , , 

e B.l  <«  + 4«) 

m 

By  differentiating  Equation  (91)  with  respect  to  r and  using  Equations 
(68)  and  (97),  the  percentage  of  heat  leaving  the  malt  without  contri- 
buting to  further  melting  may  bo  found.  The  heat  which  does  not  con- 

2 2 

tribute  to  melting  Increases  as  s /I  . For  the  SOS  diodes  used  by 
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Budeastaln,  et  al.  [2],  at  I ■ 0.56  A,  thla  amounts  to  approximately 
20Z  of  the  outflow  whan  the  melt  channel  diameter  becomes  comparable 
to  the  silicon  thickness.  For  very  short  times,  the  heat  leaving  the 
melt  channel  may  be  assumed  to  go  entirely  Into  melting  more  silicon. 
This  simplifies  the  analysis  considerably,  when  valid,  and  corresponds 
to  Goodman's  approximation  [29].  The  general  utility  of  Goodman's 
approximation  for  melt  channel  analysis  Is  limited  to  very  short  times, 
corresponding  to  those  at  which  the  physical  model  Is  least  accurate 
and  to  very  long  times  when  the  channel  Is  no  longer  cylindrical. 

As  Illustrated  In  Figure  39,  the  voltage  rises  slightly  during 
the  pulse  due  to  self-heating.  The  peak  value  may  be  taken  equal  to 
IR^,  yielding  Immediately.  If  L Is  known  from  photomicrographs  or 

electrical  measurements,  than  a can  arbitrarily  be  selected  Just 

below  the  knee  of  the  curve  In  Figure  39  and  the  corresponding  s can 

be  found  from  Equation  (95) . Substituting  Equation  (89)  Into 
Equation  (94), 


V 


1 + 


(98) 


The  value  of  6 can  be  found  by  curve  fitting.  As  shown  In  Figure  40, 
the  solution  Is  valid  for  the  falling  portion  of  the  curve  until  the 
melt  channel  diameter  becomes  comparable  to  the  thickness  of  the  silicon 
layer. 


Qualitative  agreement  with  experimental  data  by  Budensteln,  et  al. 
[2]  after  the  first  few  nanoseconds  Is  fairly  good.  While  linear 
approximations  could  have  been  used  for  X,  a,  R,  etc.  Instead  of  their 
estimated  mean  values,  the  additional  complication  did  not  seam  war- 
ranted In  view  of  the  approximations  Inherent  In  the  model  Itself,  the 
extrcsie  lack  of  good  data  on  the  thermal  properties  of  molten  silicon, 
and  the  qualitative  agreement  with  the  experlswnt.  As  pointed  out  by 
Goodman  [29],  some  Isiprovement  Is  possible  if  a cubic  solution  Is 
assumed  Instead  of  the  quadratic  solution  of  Equation  (78).  If  n 
channels  are  present,  the  total  currant  I nay  be  divided  by  n and  used 
In  the  analysis  Just  presented,  provided  that  the  separation  between 
channels  Is  large  enough  that  the  external  temperature  distributions 
given  by  Equation  (91)  have  negligible  overlap. 

Budensteln* 8 SOS  devices  are  sufficiently  thin  that  the  cylindrical 
melt  channel  quickly  reaches  both  the  surface  and  the  Interface  with 
the  sapphire,  becoming  rlbbon-llke  rather  than  cylindrical  In  less  than 
a microsecond.  In  the  following  analysis,  the  growth  of  this  ribbon 
of  molten  sUlcon  will  be  considered.  The  problems  of  large  channels 
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in  convancloMl  transistors  must  be  handled  on  an  individual  basis, 
because  such  a variety  of  possibilities  exists.  It  will  be  assumed 
that  the  solid-melt  interface  is  perpendicular  to  t’  substrate 
(Figure  41). 


For  simplicity,  constant  (average)  parameters  will  be  assumed. 
Then,  using  rectangxilar  coordinates. 


c 0 


II 
m at 


K - H(T  - T J + J*  & 
3x 


(99) 


in  which  is  the  temperature  at  the  bottom  of  the  sapphire  substrate. 

The  new  loss  term  can  be  ranoved  by  a change  of  variables,  as  discussed 
previously,  or  the  equation  can  be  solved  directly.  The  latter  will 
be  done  here. 

For  the  purpose  of  illustrating  the  origin  of  this  additional 
term,  an  average  value  is  assumed  for  the  thermal  conductivity  and 

tha  ganaration  tarn  is  temporarily  ignored.  As  shown  in  Figure  41, 
the  malt  ribbon  volume  looks  almost  like  a line  source  in  comparison 
to  the  volume  of  tha  substrata.  Than 


- Viz’  ■ Vj**  W * 
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to)  RIBBON  MBLT  CHANNEL 


to)  SUBBTRATl  AND  "UNI  tOURCE**  (INTO  MOB) 

Plgur*  41.  Depletion  region  over  eubstrete. 

By  aeklng  the  (crude)  easunptlone  thet  the  temperature  hae  reached 
ambient  temperature  at  a radius  d^^  from  the  depletion  region  and 

that  the  melt  region  can  be  approximated  by  a cylinder  of  radius  d/2, 
a rough  approximation  to  may  be  found.  Within  the  substrate, 

’•’d  - i ^ <'’d>  - 0 • 

Integrating  and  substituting  Fourier's  equation. 


36 


In  vl«w  of  Che  approximations  already  made,  an  average  thermal  conduc 
Clvlty  Is  assumed  for  Che  substrate.  For  sapphire,  It  Is  assumed 

that  ■ 0.15  tf/cm-*K.  Integrating  a second  time  and  using  the 

boundary  conditions  described  previously. 


Assuming  d - 1 urn,  d.,  - 250  um,  K - 0.15  W/cm-*K,  the  value  of  H to 

be  used  for  this  structure  Is  H > 5 » 10  W/cm  -*K.  As  Illustrated 
previously  In  Flgtire  37,  the  solid  silicon  around  the  melt  channel 
quickly  approaches  Che  melt  temperature.  It  thus  appears  that  the 
temperature  gradient  within  Che  melt  channel  always  consists  of  a change 
of  a very  few  degrees  In  a very  short  distance,  while  the  late-tlme 
taotperature  gradient  outside  the  rlbbon-llke  channel  consists  of  an 
equally  small  change  In  temperature  across  a much  wider  region.  For 
this  reason,  Goodman's  boundary  approximation  that  all  heat  leaving 
the  channel  goes  Into  melting  becomes  a fairly  good  approximation, 
unlike  the  situation  at  Intermediate  times.  Goodman's  approximation 
will  bo  used  here,  because  only  the  general  trends  are  of  Interest. 

A mors  rigorous  solution  can  be  found  by  matching  boundary  velocities, 
as  In  ths  previous  section.  The  equation  for  tamperature  outside  the 
melt  Is  replacsd  by  a far  more  compllcatad  one  containing  terms  In 
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Th«  double  aovlng  boundery  problem  erislng  from  growth  In  both 
the  poeltlve  end  negetlve  x directions  (Figure  42}  cen  be  simplified 
by  a synnetry  argument.  A single  boundary  s(t)  Is  considered,  where 
s Is  the  distance  from  the  centerline  of  the  melt  channel  to  the  melt- 
solid  Interface,  and  current  1/2  Is  used.  The  melt  channel  Is  L In 
length,  2s  In  width,  and  d In  thickness.  An  "Initial"  square  cross- 
section  of  side  a(t^)  ■ s^  > d/2  Is  asstaed  and  the  Intrinsic  channel 

resistance  la  Ignored.  Then 
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The  boundary  conditions  given  In  Equations  (68) , (69) , and  (70)  are 
still  valid,  provided  r Is  replaced  by  x.  If 


T • T - a(s  - x)  + b(s  - x)^ 
n 


(110) 


the  coefficients  are 


In  less  than  1 sec,  the  square  root  term  goes  to  one  (generation  equals 
loss);  the  solution  should  not  be  carried  beyond  this  point.  Continuing, 
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for  c.  < t < t _ . Thus,  Che  melt  chsnnsl  cesses  to  expend  eftsr 

X “ ••  uiAX 

epproxlmetely  36  msec,  et  which  elms  the  melt  chennel  hes  Increesed 
in  width  from  d Co  epproxlmetely  6.7  d.  Beesuse  the  other  psremeters 
In  Equetlon  (117)  ere  Independent  only  on  the  meterlel,  silicon,  the 
term  governing  the  mexlmum  time  is  the  effective  H,  which  wes  shown 
previously  to  depend  entirely  upon  Che  properties  of  the  substreCe 
end  the  essumptlon  of  en  Inltlel  sqiiere  cross  section  equel  to  the 
silicon  film  thickness.  The  mexlmum  width,  2s  , depends  on  the 

nAX 

totel  current  es  well.  While  the  enelysls  Is  edmictedly  crude.  It 
does  give  some  Indlcetlon  of  the  reletlve  Importence  of  verlous  persm- 

eters.  For  exemple.  If  A * 9.17  x 10^  J/kg  Is  used  (Teble  3),  C^^ 

Is  only  18  msec.  Flnelly, 


(119) 

Equetlon  (113) , with  the  epproxlmetlon  for  v given  In  Equetlon  (114) , 
cen  be  used  to  demonstrete  thee  for  I > 0.56  A,  the  entire  chennel  Is 
essenclelly  conscent  In  tempereture  et  T^.  This  egress  with  the  trend 

begun  In  Figure  36.  Figures  43  end  44  show  s(c)  and  V(c),  respectively. 
It  Is  noted  thet  Figure  40  ends  et  1 usee  end  32  V,  while  Figure  44 
begins  et  50  usee  end  50  V.  Thus,  most  of  the  ectlvlty  occurs  In  less 
Chen  1 msec,  with  e very  greduel  decreese  to  e low  voltege.  If  the  test 
pulse  Is  shorter  then  the  voltege  drops  with  the  fell  time  of  the 

test  pulse.  There  Is  reelly  no  effective  wey  to  show  the  entire  pulse 
on  e single  greph;  Figure  45  shows  e composite  of  Figures  40  end  44 
with  e log- log  scela.  A third  pleteeu  exists  et  epproxlmetely  7.7  V 
end  36  msec,  where  the  melt  channel  ceases  to  grow.  The  change  In 
growth  race  as  the  heet  loss  to  the  substrate  becomes  significant  Is 
clear ly  Indicated.  It  the  latent  heet  Is  only  half  es  large  es  Chet 
used  la  the  eeneellstleas,  the  voltege  et  the  tell  of  the  pulse  will 
effectively  be  divided  by  1.4  and  will  drop  somewhat  more  rapidly  than 
ladlceted. 


V.  cuRimiirr  amof  woono  MiAKOOiivN 

Curromt  nsde  second  braekdowa  Is  s special  form  of  breakdown 
occurring  only  In  epitaxial  planar  transistors  [4,31].  Because  this 
construction  process  Is  so  eonnon,  current  node  second  breakdown 
warrants  special  consideration.  Noreo''er,  under  certain  conditions. 

It  can  lead  to  theraal  second  breakdown  and  damage,  even  though  la 


60 


Figure  45.  Composite  voltage  curve. 


Itself  It  Is  nondestructive  [4].  Figure  46  shows  the  schematic  cross- 
section  used  by  Crutchfield  and  Motoux  to  explain  current  mode  second 
breakdown  In  an  epitaxial  planar  transistor  [32].  The  collector-base 
Jtactlon  Is  essentially  a one-sided  abrupt  junction. 

In  the  Grutchfleld-Hotoux  model  of  the  current  mode  phenomenon, 
reachthrough  is  required,  and  the  epitaxial  layer  must  be  thin.  When 
the  ohmic  drop  across  the  base-spreading  resistance  Is  sufficiently 
large,  the  emitter  will  Inject.  Avalanche  multiplication  of  carriers 
Injected  by  the  emitter-base  junction  causes  conductivity  modulation 
in  the  region  of  maximum  current  density,  with  plnch-in  enhancing  the 
effect.  According  to  their  model,  the  condition  for  Initiation  of 
current  mode  second  breakdown  is 


J - q v^,  H (120) 

e wmx  c 

In  which  J is  the  collector  current  density,  only  drift  current  is 
c 

present,  carrier  velocity  is  eeturated  at  a maximum  value,  and  the 
majority  carrier  density  la  equal  to  the  doping  density  of  the 

epitaxial  layer.  Crutchfield  and  Motoux  then  use  Fletcher's  [33] 
solution  for  the  current  density  at  a distance  r from  the  center  of 
the  emitter: 
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Figure  46.  Schemeclc  of  eplcaxlal  planer 
transistor,  cross-sactlonal  view. 
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where 


tg  ■ radius  of  emitter 
Jg^  ■ current  density  at  tg 

* large-signal  coonon  base  current  gain 
Og  • base  conductivity 
Wg  • base  width 

Jg  ■ average  currant  density  In  emitter 

end  integraee  to  gee  enltter  current  I_  required  to  Initiate 
■ode  second  breakdown:  ^ 

Je 


(121) 


(122) 


current 

(123) 
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Th«  preceding  equeClone  are  repeated  here  because  typographical 
errors  appear  In  the  original.  The  authors  claim  that  a sudden  change 
In  conductivity  modulation,  l.e.,  the  current  mode  second  breakdown, 
takes  place  at  the  critical  current  defined  In  Equation  (123),  but 
they  do  not  say  why. 

Rabum  [4]  provides  a more  meaningful  explanation  of  current  mode 
second  breakdown  by  an  amplification  and  clarification  of  the  Ideas 
expressed  In  the  Grutchfield-Motoux  model.  He  also  points  out  that 
their  limitation  to  thin  epitaxial  layers  Is  not  necessary  and  lists 
criteria  for  determining  whether  current  mode  second  breakdown  will 
take  place  in  a given  structure.  His  analysis  begins  with  the  same 
differential  equations  as  those  found  In  Section  II,  but  he  does  not 
neglect  mobile  carriers  In  Poisson's  equation  (shown  here  for  cartesian 
coordinates) : 


in  which  the  mobile  carriers  sre  assumed  to  be  flowing  only  by  drift 
and  at  their  saturation  velocities: 


■J,  ■ 1 " ’tp 


J • q n V 
n 


In 


An  Injection  ratio. 


J (W) 


(125) 

(126) 


(127) 


Is  defined  and  results  explained  In  terms  of  It.  Figures  47,  48,  and 
49  show  electric  field  variation  serose  the  depletion  region  for  a 
one-sided  abrupt  collector-base  Junction  at  reachthrough  and  y • 0.1, 
0.5,  and  1 [4].  The  distribution  Is  seen  to  vary  greatly  at  higher 
values  of  J. 

The  following  la  e brief  sketch  of  the  Rebum  model.  An  emitter 
configuration  at  reachthrough  Is  considered.  As  current  through  the 
base  Increases,  the  voltage  across  the  base  spreading  resistance  will 
eventually  reach  the  avalanche  breakdown  voltage  of  the  emitter-base 
Junction.  This  will  occur  first  near  the  top  surface  of  the  chip. 
Once  current  flows  through  the  avalanching  region.  It  can  proceed 
down-ward  to  the  planar  central  portion  of  the  emitter-base  Junction 
depletion  region  (which  will  be  forward  biased),  where  It  will  be 
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Injsctsd  into  ths  bass,  carrisd  across  tha  base  with  soma  aid  by  the 
build-in  field,  and  go  onward  to  the  avalanching  collector-base 
Junction,  as  proposed  by  Crutchfield  and  Motoux.  The  Increased 
injection  caused  by  the  subsurface  avalanche  and  consequent  sudden 
change  in  injection  ratio  are  the  missing  points  in  the  Grutchfield- 
Motoux  model.  As  pointed  out  by  Rabum  [4],  the  avalanching  subsur- 
face region  may  carry  so  much  current  that  it  goes  into  thermal  second 
breakdown  and  is  damaged.  Doping  densities  are  high  enough  that 
tunneling  may  be  present.  According  to  Rabum* s model,  current  mode 
second  breakdown  will  occur  when  total  currant  - \l'^  in  the 

configuration  and  when  I2  ■ ^br^’^b  ^ ^cbo  vhere 

is  the  built-in  voltage,  and  r|^  is  the  base  spreading  resistance. 
Rabum  checked  the  model  by  eliminating  r|^  from  the  two  equations  and 
comparing  experimental  and  predicted  values  of  I2: 


(128) 


Because  ■ 20  1^,  the  velue  of  could  not  be  cucpected  to  remain 

completely  constant;  nevertheless • the  results  were  In  reasonably 
good  agreement. 

It  may  be  possible  to  extend  the  analysis  somewhat  by  replacing 
the  base-spreading  resistance  with  a Junction  field  effect  transistor 
(JFET)  as  Illustrated  In  Figure  50.  Diode  D^  represents  the  avalanching 

subsurface  region,  and  diode  represents  the  collector  In  reach- 

through,  so  avalanche  diodes  ("Zeners")  are  used  In  the  equivalent 
circuit.  Diode  D_  Is  an  ordinary  forward-biased  diode,  and  R Is  the 

A S8 

substrate  resistance.  Each  symbol  represents  a nonideal,  lumped  element; 
resistances,  capacitances,  etc.  are  not  shown  explicitly.  Standard 
JFET  analysis  may  yield  a much  better  expression  for  the  nonlinear  base 
spreading  resistance.  It  should  be  noted  that  Figure  45  Is  somewhat 
misleading  In  that  It  does  not  warn  that  a ring-dot  geometry  Is  Involved 
and  that  the  effective  source-drain  current  flows  radially  In  a 
cylinder. 
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Figure  SO.  Equivalent  circuit  with  nonideal  lumped  elements  (pnp). 
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